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ABSTRACT 


The development of a nonplanar lifting surface method having a 
continuous distribution of singularities and satisfying the tangent flow 
boundary condition on the mean camber surface is given in this 
dissertation. The method predicts some incompressible longitudinal 
aerodynamic coefficients of rectangular wings which have circular-arc 
camber. The solution method is of the integral-equation type and the 
resulting surface integrals are evaluated by either using numerical or 
analytical techniques, as are appropriate. 

Applications of this method are made and the results compared 
with those from an exact two-dimensional circular-arc camber solution, 
a three-dimensional flat-wing solution which represents the camber by 
a projected slope onto the flat surface, and a flat-wing experiment. 

From these comparisons, the present method is found to predict well the 
flat-wing experiment and limiting values , in addition to the center 
of pressure variation at an angle of attack of zero for any camber. For 
wings having camber ratios larger than about 1.25# and moderate to high 
aspect ratios, the results of the present method deterioriate due to the 
inadequacy of lifting pressure modes employed. 
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VI. INTRODUCTION 


Solutions for the steady aerodynamic characteristics of wings 
began with Prandtl (ref. l) and his famous "lifting line" concept. In 
this concept, the loading was concentrated along a straight line and 
all induced effects on the wing came from the trailing vortex sheet. 

The local effective angles of attack were composed of the geometric 
angles and the induced angles along the lifting line. This concept 
was useful in predicting the characteristics of lifting systems 
employing rectangular wings of high aspect ratio by representing the 
span loading as either an averaged constant value or a semi -elliptical 
variation. However, as the planforms became more complex, additional 
improvements in the theory were necessary in order to predict the 
aerodynamic characteristics. The improvements can be classified in 
three ways: (l) extension "of Prandtl*s concept to other planforms, 

(2) inclusion of the induced-camber effects in the lift distribution, 
and (3) allowance for nonplanar lifting systems (other than biplanar). 

Prandtl 1 s "lifting line" method and the extensions employed to 
allow for the consideration of other than rectangular planforms provided 
span load distribution , and in some instances , determined the induced 
drag as well. For instance, Multhopp (ref. 2) devised a method for deter- 
mining the span load distribution on any< wing that could be represented 
by a straight lifting line. Mutterperl (ref. 3) developed a method for 
determining span loading on sweptback wings by using a bound vortex 
(that is, a lifting line of variable strength) lying along the quarter 
chord of each wing panel, as well as a sheet of trailing vorticity. 


1 
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The strengths of the vortex system were constrained to give tangent 
flow on the wing three-quarter chord line. This relationship between 
the bound-vortex position and the tangent-flow-boundary line is the 
same as for the two-dimensional solution, which is one that develops 
the correct lift value. Weissinger (ref. b) developed both lifting- 
surface and lifting-line methods for sweptback wings, and used the 
tangent flow constraint at the three-quarter-chord location to obtain 
solutions. 

The inclusion of the induced -camber effects on three-dimen- 
sional planar wings is seen in two important ways. One is the removal 
of the constraint that the section center of pressure be at the 
quarter chord, and the other is the improvement in predicting the 
lifting-pressure (difference between lower and upper surface 
pressures) distribution because of the increased amount of information. 
All of this is accomplished by either expanding the number of chordwise 
mode shapes or placing additional bound-vortex elements and control 
points along the chord, or, as Faulkner (ref. 5) did, with a combina- 
tion of the two. The lifting surface method of Faulkner’s is called 
"vortex lattice theory" and is based on representing the induced 
effects of the lifting system by using a lattice of horseshoe vortices, 
the strengths of which are determined by a series of modal functions. 

In this solution, there would be as many coefficients in the series as 
control points (points at which the tangent flow boundary conditions 
are satisfied), but not as many as there are horseshoe vortices. Also, 
Multhopp (ref. 6) developed a lifting surface method called "subsonic 
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lifting surface theory" which used two modes of pressure chordwise and 
two control points at several spanwise stations and could be applied to 
arbitrary planforms. 

Since reference 6, there have been many other lifting surface 
methods published for steady subsonic flow. Each new method, in 
general, begins with an already-developed fundamental method and seeks' 
to improve its accuracy by (l) adding additional chordwise loading 
terms (refs. 7 and 8) , (2) reformulating the solution (ref. 9 ) »' 

(3) performing the integrations differently (ref. 10)., or (4) extending 
the original concepts. For example, extensions to Faulkner's work 
can be made easily by requiring that the distribution of vortex 
strengths be individually determined (not expressed in a series) by a 
solution to a set of linear simultaneous equations employing as many 
control points as unknown vortex strengths. (See, for example, 
references 11, 12, and 13.) Many of the latter-re ferenced methods 
would not have been practical to use before the advent of the high- 
speed digital computer. 

Solutions for nonplanar lifting surfaces, such as wings with 
dihedral, have been made with the vortex- lattice method (refs. 11 and 
lU) , kernel-function-integral equation (ref. 15), and an asymptotic- 
expansion procedure (ref. 16 ). Other nonplanar solutions are avail- 
able for intersecting tail surfaces j among them references 11, 17, 
and 18. In each solution, the boundary conditions are satisfied on 
the chordal plane. However, for wings which are nonplanar in the 
chordwise direction (that is, for a cambered wing) , simplifying 
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assumptions are usually made in order to effect a solution. The major 
assumption is that of satisfying the boundary conditions of the cambered 
wing on a lifting surface lying in a plane formed by the longitudinal 
and spanwise variables of a wind axis system by projecting the wing 
slopes onto the plane. (See sketch 1 and, for example solutions, see 
references 8 and Ik . ) 



By doing this , the effect of the vertical displacements between 
influencing and influenced points is not taken into account. Another 
method (ref. 19 ) , which satisfies the boundary conditions on the surface 
of the cambered wing, does not satisfy the Kutta condition at the 
trailing edge because the solution is based upon solving a series of 
two-dimensional cross-flow problems. The formulation of the vortex- 
lattice method, as given in reference 11, appears to be general enough 
to account for vertical displacements as well as to satisfy the boundary 
conditions on the surface. However, constraints would have to be 
imposed in order to guarantee that the trailing vortex system only 
exited the wing at the trailing edge. 

In consideration of the limitations of the methods just dis- 
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cussed, it seems appropriate to develop a solution for the cambered 
wing problem which would (l) satisfy the boundary conditions on the 
surface, (2) account for the vertical displacements between influencing 
and influenced points, and (3) satisfy the Kutta condition at the 
trailing edge. This dissertation describes such a solution for 
rectangular wings having circular-arc camber (see sketch 2). 

(f) 


A 

v,x 

1- Z 

Sketch 2. Representation of a rectangular wing 
having a constant circular-arc camber. 

The solution developed herein uses a more general form of the subsonic 
integral equation than the methods discussed for planar wings, because 
it must relate the normal component of perturbation velocity - rather 
than the downwash - to the lifting pressure, and it must include the 
effects of vertical displacement. Furthermore, the present method of 
solution employs both numerical and analytical techniques. The latter 
are applied in a small region (square box in the surface) surrounding 
the control point, because it is at the control point that the integral 
equation has a singularity . This technique encompasses expansions in 
small parameters, followed by analytic integrations over the box. 
Outside of the box, the integral equation is evaluated numerically. 
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Results obtained with the present method are compared with both 
an exact two-dimensional circular-arc-camber solution and a three- 
dimensional flat lifting surface which employs the slope-projection 
procedure to account for camber. Due to the general unavailability of 
experimental data for circular-arc-cambered wings , comparisons with 
experiment are only made for flat wings. From the comparisons made, 
conclusions are drawn about the applicability of the present method 
over the camber- and aspect-ratio and angle of attack ranges. In 
addition to these applications, the present method may be extended to 
certain wing-flap systems that approximate the circular-arc-camber 
constraint. 

Three appendices are given: appendix A discusses some of the 

integral types which must be integrated over the box surrounding the 
control point; appendix B presents the equations used to compute some 
section and wing aerodynamic coefficients; and appendix C gives the 
development of the aerodynamic characteristics of the two-dimensional 
circular-arc airfoil, included herein for completeness. 



VII. THEORETICAL DEVELOPMENT 


Basic Formulation 

Concepts 

The method of solution for pressure loading over a wing 
having circular-arc mean camber line (sketch 2) will begin with 
Poisson's equation (see ref. 20), as does the flat-wing solution. 

The basic concept here is that the general perturbation velocity 
field will be the same for both wings, with differences occurring 
primarily in the vicinity of and aft of the wing. It is possible, 
therefore, to make use of some of the solution techniques which have 
been developed for the flat wing in seeking a solution to this 
problem. 

From reference 20, the three-dimensional Poisson's equation 
can be written for the perturbation velocity in the freestream direc- 
tion as 


uU, y ,z) =-h///?v 2 udv ♦ fer //[u-|j (i) - i fa] 

\T <3 ' / J 


dS 


( 1 ) 


By representing the wing surface with a sheet of pressure doublets, the 
perturbation velocity can be related to a velocity potential. The 
velocity potential due to each pressure doublet is of the form 
-( z-z-i) 

3 — and is a solution of Laplace's equation 


V 2 <(>(x,y,z) = 0 


( 2 ) 


7 
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Upon relating u(x,y,z) to <(> by 

u(x,y,z) = <t> x (3) 

substituting equation (3) into equation (2), and interchanging orders 
or differentiation, it follows that 

y 2 u = 0 (4) 

Hence , 

u(*,y,z) = -fc //[ » ■£ (i) - i £] as (5) 

The surface S includes both the upper and lower surfaces of 
the wing, as well as the wake upper and lower surface area. Since 
the flow is assumed to remain attached to the wing, only leaving at 
the trailing edge and thereby satisfying the TEutta condition, the 
perturbation velocity in the freestream direction on either side of 
the wake surface must be the same. Consequently, the lifting pressure 
is zero. Hence, upon examining the wake portion of equation (5)j 



it can be seen that , since ^ ^ and on the upper surface 

is opposite to the same tern on the lower surface, the first terms in 
each integral cancel. That is, 

//K - - J I? (r)] 158 = ° 

S Wu 


( 6 ) 
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Also, for wake shapes that are straight or circular arc, the normal 
derivatives of u are continuous across the wake because they are 
related to the continuous streamwise derivatives of the normal velocity 
components through irrotationality . This leads to cancellation of the 
second terms in each integral. That is, 


* //[-7^ ] dS = 0 


( 7 ) 


w 


u 


*1 


Thus , the integral equation that remains to be evaluated is 

s u s i (8) 

The coordinate system chosen in which to solve this equation 
is cylindrical polar because of the ease in applying the boundary 
conditions* In this system* the mean camber surface lies on a surface 
of constant radius. Hence, the normal direction is along a radius 
vector R*, the tangential direction is along $ , and the spanwise 
direction is along y^. This leads to 



u(x ’ y ’ z)= c/Im -(i)| (3* - y 


R*=R 


R*=R 


R*=R 


dS 


( 9 ) 


From sketch 3, 
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Sketch 3. Velocity diagram on the surface 
of a circular-arc section at R,t? 


it can be seen that 


% = % sin * “ % cos * 


(10a) 


and 


u x = u^ sintf - u n cos d 


(10b) 


By differentiating equations (10a) and (10b) with respect to R*, the 
following equations are obtained: 


3u 


3u 


u 


= sm 




U u 


- cos 




3u 


n 


3R* 


3R* 


3R* 


(Ha) 
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and 


9u 

9R* 


1 = sin' 3 


9 “t .85: 

• 1 - cos 3 — — 

9R* 9R* 


(lib) 


with u n being continuous across the plane of the wing, or along 
R* . 

The irrotationality condition, when written in cylindrical- 
polar coordinates , gives the needed relationships between the per- 
turbation velocities tangent and normal to the camber surface as 


9u tu _ _ 1_ 9u n _ ut u 

9R * R* 9T ~r* (12a) 

and 

3u t x = 

9R* ~R* dd R* (I2h) 

Therefore, when the results of equations (lla) , (lib), (12a) , and (12b) 
are substituted into equation (9) , the result is 


u(x,y ,z) 


1 



(u_(. -U|.^)sini? 



Let 


- sin 0 



+ *t 
rR* 


i 


> 

') 



|r*=r 


(13) 


Y c = (ut u - u ti ) 


(l^a) 
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Then 


u(x,y , z) = 




(lfcb) 


The velocity potential associated with this perturbation velocity can 
be determined by 


4>(x,y,z) = f* i u(x' ,y,z) dx’ (15) 

-CO 

The minus infinity value is chosen as the lover limit because there 
the velocity potential is zero. 



Sketch 4. Path of first level of 
integration. 


From sketch 4 it can be seen that 

x ! = R* cos0* (l6a) 

and hence 


dx' = dR ! cos0 ! + R ! (“Sin0'd0 f ) 


( 16 b) 
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The x' integration must take place at a constant z elevation, 
as in the flat-wing solution; therefore, by setting 

z = R' sin0' = R sin0 (17) 

equation (l6b) can be re-expressed in terms of the constant z as 

dx* =d(z/sin0') cos0 ' - zd0' (l8) 


Since 


dR' = 


dR' / JQ ,V_ d(z/ sin0 ' ) 

d0 * \ aw / _ d© 7 



-ZCOS0 1 

sin 2 0 ' 



(19) 


then 


dx f 


-zcos 2 0 1 
sinks' 



zd0' 


-zd0 ' 

sin^e ' 


( 20 ) 


or 

j t _ -R 1 d0 * -Rsin0d0 ' 

sin0 ' sin 2 0 ' ( 21 ) 


and the limits of integration go from tt to 0. Therefore, 


<J>(R,0,y) 



-Rd0 * 
sin 2 0 ' 


( 22 ) 


Once the velocity potential is known, differentiation with respect to 
R will yield the perturbation normal velocity. Hence 


"n 


_ 3(t>(R.0.v.). 

3R 


isin0 


(R.0.y) 


(23a) 



Ik 


and so 


<= fl D> I 


^ (R,0,y) = 


-Y c siir8Rsin0 f 5 / 2 _ 


sin^e' 


W (r) + (FF)| I d0'Rd ^dy^ 

R*=R R*=Rj/ 


(23b) 


After performing the indicated differentiation, equation (23b) can be 
written as 


(HjQjy) 




Ry e sin^sin 0 f _3^_ /j\ /1 \ 9 / 1 \| 1 ) 

sin 2 © ' 3R3R* \r/ \R*/ * ( r / / d6 ' 

L R*=R R*=R |r*=r -I / 


1 Rd$ dy x 

(24) 


or 


*Ifc 


u n (H,9,y) = f-l /-Ry c si^sin0 


,..n 


j.i 


sin 0' 


[■55* (r) * (ihi 

L R*=R 


R*=R 


R3 2 

3R9R* 


(*. 


r) + ® (I? lj de ' dS ^ 

■ R*=R I R*=R J ) 


(25) 


The total normal velocity of a point lying on the wing surface ( R =R) 
must be zero in order for the flow to be tangent to the surface; 
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hence , 


lim u n (R,0,y)' + V cos0 = 0 (26) 

R+R 


Since 

V = Ap . 2V _ AC P V (27) 

c pVsint? 2V 2sint? 


equation (25) can he re-expressed as 


lim "n 
R-*R Ysin0 


= -cot0 


_ lim 1 f f 

R->R Tnr J J 


RAC, 




Now r (see sketch 5) is defined as the straight line distance 
between a point on the pressure doublet sheet at (R*,0, y^) or (x^,y^,z^) 
and a point on the path of integration (R',0',y) or (x' ,y,z) . 

r = •^(x 1 - x' ) 2 + (y~ -y) 2 + (z x - z) 2 (29a) 

where 

x^ = R*cost> (29b) 


x'= R'cosO' = R sinBcoxS ' 


(29c) 
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H 

II 

H 

(29d) 

y = y 

(29e) 

= Resin'S 

(29f) 

sin0 f = Rsin0 

(29g) 



Sketch 5. Distance between point on pressure doublet 
sheet and a point on the path of integration 

For R* = R, the sheet of pressure doublets is coincident with the wing 
surface, and for R-tR, the field point moves to the surface since this 
is where the boundary condition of tangent flow must be satisfied. By 
performing the indicated partial differentiation and then taking the 
limit as R ->R, equation (28) becomes after simplification 


IT 


cot0 


lim R_ (( 

■ y " R-HR 8irJJ ' 

*0 


(y 1 - y) : 


- 2zsinS 


J T\ 


-2zcos$ I SSS^ASl * 3z 


r 3 sin 3 0' 


(/: 


d0 1 


r 5 sinV 


Hi FFPF 


(Rzsin 2 S) 


J £OS0^d0_[_ f _2Rz sin 1 ® co + (y^ _ y) 2 cossj 

JL r 5 sin 3 0 ’ *■ 

h / r 5 sin 2 0’ [ Rz(cos2,9 " sin ^) + (yi - y) 2 sins]^ 


dSdy 1 (30) 


First Level of Integration 

Each of the 0' integrals in equation ( 30 ) can he integrated 
directly by employing the transformation 

p = cot0' (31) 


which results in 


H 

ti 

cr> 

CD 

RcosS - zcot0 

J r 3 sin 2 0 1 ztCRsin^ - z) 2 + ( y ^ - y) 2 ] 

/D + Bcot0 + Ccot 2 0 


(32a) 



cos6 T d6 f 
r 3 sin 3 0 ' 


^1 

z 2 [(RsinS - z) 2 + (y^ 


tft* 


iRzcos'ScotG - R 2 




+ 2RzsinS - z 2 - (yq - y)‘ 
+ Bcot0 + Ccot 2 0 


+ RcosS 


(32b) 
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d0 1 


r 

J r 5 sin 2 0' 3z[ (Rsind-z) 2 + (y-^-y) 2 ] 

J ^ + ■) - 

\D+Bcot0+Ccot 2 0 (Rsin'fl-z) 2 + (yq-y) / 


Reps'# -zcot0 
A)+Bc ot 0+ C c o t 2 0 


(Rsin#-z) 2 + (y 1 -y) : 


/* 0 /*( 

j cos9 T d0 * _ 1 + Rcos# j 

J r 5 sin 3 0 ' 3z z (D+BcotO+Ccot*©) 3 /2 z J 


d0 ' 

r 5 sin z 0 ' 


and 


/: 


d0 1 


r 5 sin 4 0 ' 


^r^r + H- E!+ "- (y i- y)2 ] f 

z J r sin 0 ’ z L J J 

it 77 

jcos'S | 

J 


d0 1 


5 . 5 

r sin 


+ 2Rzcos^ j cos0 T d0 * 

~ z 7 / r 5 sin 3 0 f 
t r 


where 


= ^ + Bcot0 f + Ccot 2 0' 


and 


D = R 2 - 2Rzsin# + (y - y) 2 + z 2 


B = -2Rzcos , 9 


C = 


z 


2 


With the substitution of equation (32) into equation (30) and the 
resultant simplification of terms, the resulting equation is 


(32c) 

(32d) 


0' 

(32e) 

(33a) 

(33b) 

(33c) 

(33d) 
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totei -Kfi-rr M 

1 y R-*R 8 tt J J p 

S 

. i 

[(y-i-y) 2 1 # 

L R 

R(cos^-cos0) n 

z[ (Rsind-z) 2 H 
2cos^ 

► ( yi -y) 2 I | 

|r 2 cos 8 cos0- 

/D+Bcot0+Ccot z 0 

- > R^-R 2 sin 2 0+2R 2 sin0sin'3-(y2-y) 2 

z[ (Rsin’fl-z) 2 + (y 1 ~y) 2 ] 

[ v*l)+Bcot0+Ccot z 0 


+Rcosd 


+ |2R 2 ^sinSCcos 2 ^ -sin 2, 9 ) + sin (sin 2 0-cos 2 ^ )) + (y-^-y) 2 * 


> 

cos^Xl 

r 

R(cosd-cos0) 

\ sin0 /J 

✓ 

JJ (Rsind-z ) 2 + (y x ~y ) 2 ] 

^)+Bcot0 +Ccot z 0 


( l 

'D+Beotfl+f 


+ (y x -y) 2 ) 


VD+Bcot0+Ccot 2 0 (Rsind-z) 2 + (y n -y ) 2 / (Rsind-z) 2 
+ |2R 2 sindcos (sin^-sin©) + (y-^-y) 2 cos-i‘|* 

d^dy-L 


+ (yi-y) 2 ]]^ 


z(D+Bcot0+Ccot 2 0) 3 /2 


(3U) 


Pressure Functions 

The unknown AC p distribution appearing in the surface integral 
of equation (3*0 is prescribed herein (see sketch 6) to be composed of 
up to five terms in the Birnbaum series chordwise with the undetermined 
coefficients being functions of spanwise position. 



Spanwise mode shapes 


VXAW 

0 7T/2 i r 

€ 

Chord wise mode shapes 

Sketch 6. Loading mode shapes. 


ac p'^i> 


N-l 2 


EE 

n=0 m=0 



even 


Sketch 7. Relationship between e and * 
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Final Form of Integral Equation 

Since the chordwise pressure modes are defined in terms of a 
variable, e , which goes from zero at the leading edge to ir at the 
trailing edge, it is necessary to change the chordwise variable used 
in the surface integration from *9 to £ . Sketch 7 shows the 
equivalent edge locations in the e and *9 coordinate variables. 

The general expression relating £ and *9 is 


and so 


£ = cos 


-1 


1 - 


2R 

c 


cos(*9 + a ) - cos(0 


o + “ >] 


(36a) 


d a 


c sin £ de 

2R sin(-d + a ) 


(36b) 


Hence, by combining equations (3*0, (35a), and (36b), the following 
equation results : 


cot0 


lim 

R-*R 


n=0 m=0 j,/2 \ ** 

even u/c: 


0 f 

j ^ sinededy^ 


sin(a+a) 

(37a) 


which can be written in a more convenient form as simply 


cot0 


lim c 


c f n m «/I^ r [ h (e)(| sln ^f dy i 

'TT j / 4 11111 1 I n I I sin (‘8 +a) 

n=0 m=0 J J w 

even -b/2 0 


y r->r T&n 


^See equation (3*0 for this expression. 


(37b) 
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Since the 9 f integrals have been determined and the pressure 
mode shapes specified, the surface integration can now be considered. 
This is accomplished, in general, by means of numerical quadrature. 

The quadrature basically employed is Gaussian in both directions with 
the exceptions (as shown in Sketch 8) over the spanwise range which 


Number of Integration Steps 
Region Chordwise Spanwise 

I 

I r 50 Gauss 30 Gauss 

M 100 Gauss 70 Gauss 

U 200 Gauss 51 Trapezoidal 

3E 




Sketch 8. Integration regions and number of integration 

points used in each for low-to-moderate camber. 



2h 


correspond to the width of the square boxes centered around the control 
and its image point, the image point only occurring when 

e < i8o° - e < 0 . . 

O “ "1 

In this spanwise range , 51 evenly spaced points are employed in the trape- 
zoidal rule integration, and Gaussian quadrature, which employs 200 
steps divided according to arc length ahead of and behind the boxes , is 
used chordwise. The number of Gaussian steps and trapezoidal rule points 
have been established by convergence studies for an aspect ratio 20 
rectangular wing with a box length of .01 chord and a camber ratio of 
.00125 at a = 5°. This set of integration steps and box length have 
also been found to be adequate for cambered wings of other aspect ratio. 

Because numerical integration cannot be carried out across the 
high-ordered singularity which exists at the control point, of the form 

Az 2 5 this area is split off from the general surface 

integration and is treated separately. The image point, which is the 

cambered wing equivalent of the flat wing singular strip, has a 

singularity of the form 1 and [ 1 | 2 ; hence, the 

Az z + Ay* [ Az* + Ay Y ] 

box surrounding it must also be broken off and treated separately. The 
procedures employed are described in the following sections. 

Surface Integration Over the Control-Point Box 

The surface integration over the box is accomplished by: 

(l. establishing basic expansion variables for small spanwise and 
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angular displacements away from the control point; (2. using these 
expansions to approximate each fundamental form which appears inside 
the braces in equation (3*0; and (3. using finite part concepts, 
integrate the resulting combination of forms over the box. 

Expansions 

The expansion variables in which the expansion process is to be 
conducted are 

y 1 = y + to (38) 

and 

* = 0 + a ■ (39) 

where u) and a are the small positive displacements in the spanwise 
and angular directions, respectively. The angular displacement is 
chosen in terms of rather than e (the coordinate in which the 
loading is specified) in order to simplify the expansion process. 

Approximate Forms 

By employing the above expansion variables, the following 
approximate forms are obtained near the control point for the general 
expressions given in equation (3*0. They are to 0(w 2 ) and O(o 2 ): 

( y l~ y ^ - 2Rsin9sind + 3Rsin 2 d = ~ + R(sin 2 0 - 2a 2 sin 2 0 + 

R R 


3cf 2 cos 2 0 + Uasin0cos0) 


(UOa) 
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(Rsinfl - z) 2 + (y^-y) 2 = R 2 O 2 cos 2 0 + to 2 (HOb) 

D + B cot0 + C cot 2 0 = R 2 a 2 + to 2 (l*Oc) 

2 

R(cos$ - cos0) = - R cos0 + asin© (UOd) 

cos'-d = ^1 - ^ cos0 - Osin0 ( i+Oe ) 

R 2 cos , 3 cos0 - R 2 - R 2 sin 2 0 + 2R z sin0sind - (y 1 ~y) 2 - 

- [R(a 2 /2)(l+sin 2 0) + a) 2 - R 2 asin0cos0] (Hof) 


2R 2 [sin0(cos 2, 3 -sin 2, 9 ) + sin8 (cos 2 0 - cos 2 '3 )] = 


(-0COS0 + 0 2 sin0 - 0 2 sin 3 0) 2R 2 


(Hog) 


(72-7) 2 sin^ -sin0 


+ cos 2 ^ 
sin0 


w 


cos 2 0 

sin0 


(UOh) 


2R 2 sin 1 8 cosfl (sind - sin0) + (y 1 -y) 2 cos’8 = 2R 2 [0sin0cos 2 0 

+ 0 2 ( cos 0 - sin 2 0cos0)] + U) 2 cos0 (Uoi) 

These approximate forms are generally applicable for all angular locations 
of control points between 0 O and 0^ . An exception to this is at 
0 = 90° where the approximate form of (RsiniS - z) 2 + (y^ - y) 2 , that is, 
R 2 0 2 cos 2 0 + w 2 , tends to simply w 2 . If more terms are retained in the 
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original expansion of the above expression a more accurate result is 
obtained and it is 

R 2 |J- sin 2 0 + cr 2 cos 2 0 - O 3 sin0cos0 + w 2 

It can be shown that this more general form tends toward the original 
form for - cotdj 2 « 2cot 2 0 and tends to R 2 ^-— + w 2 for 0 90° , 

a much different result than obtained previously. 

Both approximate forms, when used, would appear in the 
denominator and hence be singular in a . However, the more general 
form would be singular in a to an order higher than that encountered 
with the original form (a 4 vs. a 2 ) and hence lead to new expressions 
to be integrated. These integrations show that the box contribution to 
be just the 0 = 90° limit of the 0 not near 90° solution (original 
approximate form) whose development follows. 

From a practical standpoint 0 = 90° should be avoided as a 
control point location so that numerical difficulties outside the box 
associated with the higher ordered singularities can be circumvented. 

It has been found that in order to assure convergence of the numerical 
integrations j© - 90°j should be less than 1.5° • 

Integrated Results 

With these approximated forms substituted into equation (34), 
the surface integration over the box becomes 
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l_im R 
R+R 8tt 


0) CJ 

f f Ac p (0 ’y) \f 


+ R(sin 2 0 - 2a 2 sin 2 0 + 3a 2 cos 2 0 


-u) -a 
-i UasinBcosoJ 


- 1 P R[ (a 2 /2)cos0 + CTsin0] 1 

;(R z O<cos 2 0 + a,*) [ +1 J 


[ R 2 (a 2 /2 ) (l + sin 2 0) - R 2 asin0cos0 + (i) 2 
/ R 2 a 2 + w 2 ' 


-(2 - a 2 )cos0 + 2osin0 
z(R 2 a 2 cos 2 0 + a) 2 ) 


+ R((l-a 2 /2)cos0 - Osin0^j 


+ [ (-acos0 + ^ a 2 sin0 - a 2 sin 3 0)2R 2 + a) 2 cos 2 0/sin0] • 


1 |r[ ( a 2 /2)cos9 + Osin0] , 

ni l . 2 \ 

[RVcos^e + L y RV + - u ,- 1 

^R 2 a 2 + co 2 R^a^cos^O+a/ J 


R 2 a 2 eos 2 0 


+ UTCOS0 


+ j + ^ as i n 6 cos2 9 + CJ 2 (cos0- 5/2 sin 2 0cos0)j 

1 __ 2 1 (hi) 

J z(R0 2 + w 2 ) 3 / 2 j 


where 


RO = 0) 


(42) 


This integral expression can be integrated with any standard set of 
integral tables, such as reference 22, and the result is to O(o°) (to 

be read order of a to the zeroth power, i.e., a independent terms) 

fSee appendix A for an example of the manner in which the 
integrations must be conducted. 
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i+sinG £l - £2|ii + £2|!£ * * * J AC p (0,y) (1*3) 

where 

(b/2) V 1 - (b/2) h n (e 0 ) ihk) 

There are some interesting features of this box integration result. The 
first is that the lowest ordered non-zero result does not depend on the 
box size. The second is that, like the flat-wing solution (unpublished 
analysis performed by Mr. P. J. Bobbitt), the leading term is determined 
to be basically a ^AC^(0,y) . Thirdly, the cosine terms can be 
thought of as an 0(7°) camber correction to the flat-wing solution 
(0 = 90°), which becomes more important as the camber increases. 
Fourthly, the O(^) terms which are omitted can also become important 
for 0 angles near 0° and l80° because there the o(^°) terms 
approach zero. However, since the effect of the 0(^°) terms can be 
controlled by selecting small values of 0 , and since some of the 

O(o) terms were obtained by expanding about 0 = 90° , the results 

presented herein are only accurate to 0(^°) • Hence, solutions for 
very large camber ratios, which lead to 0 angles near 0° and l80° , 
are attempted with caution - as should those employing the slope- 
projection technique. 

Surface Integration Over the Image-Point Box 
The same procedures used in arriving at the expanded form of the 


N-l 


v e ’ y) = L L 


n=0 m=0 

even 
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general expression are employed about the image point, with the exception 
that, instead of expanding about 9 , the angular expansions take place 
about l8o° - 0 . This only changes the signs on the cos0 terms 
which result from the terms sin(l80° - 0 + a) and cos(l80° - 0 + c). 
When these changes are followed through the simplified integral equation, 
integrands which yield the O (o°) terms are determined to be identical 
to those for the control point box, but they do have slightly different 
coefficients. Thus, the image point contribution to the normal induced 
velocity at the control point is to O(o°) 



AC p (l 80 °-e,y) (U5) 


This is an interesting result in that for control points in the first 
quadrant (image point ahead of control point) the term ^1 + 
becomes 2 , and for control points in the second quadrant (image point 
behind control point) the term (l + j qqsQ ~| 

of this term is identical to that of the modified kernel function in 


/ n , cos0 \ 

(1 + ”1 qt 1 is zero. 

\ COS0 J 


The behavior 


the flat-wing solution at the spanwise location of the singularity. 


Solution Technique 

The solution of the subsonic integral equation given in 
equation (3^) is made by requiring that the flow be tangent to the 
surface at a number of control points. There are as many control 
points as there are pressure mode-shape combinations. For each control 
point, the surface integrals are evaluated either numerically or 
analytically for each mode-shape combination. The sum of these 



integrations is stored along a row of matrix A. The unknown coefficients 
of the pressure mode-shape combinations, c^ , are determined after 
both the boundary condition value and influence coefficient matricies 
are filled by solving the matrix equation 



VIII. RESULTS AND DISCUSSION 


Gener al 

In this chapter the results of calculations made using the 
present method with various box sizes , control point locations 
and number of control points are analyzed followed by a discussion of 
results for specific wings. The first wing results presented are for 
a flat rectangular planform for which experimental data are available 
and are of interest primarily for the check it provides of the present 
method. Cambered wing comparisons are only made on a theoretical basis 
because an extensive literature search uncovered only one paper (ref. 

22) containing experimental data for cambered wings of the type 
considered herein. Unfortunately most of these data were obtained with 
the wing having some regions of separated flow. This happened because 
the airfoil shapes used were developed from thin plates which had only 
simple rounding at the leading edge and hence did not promote flow 
attachment. The present method is developed for attached flow conditions 
and hence not applicable to these data. 

The theoretical results used to compare with those of the 
present method come from methods which account for camber by the slope- 
projection technique with the lifting surface constrained to the X-Y 
plane; in particular, the methods of Scholz (ref. 23) and Multhopp (ref. 
8). Comparisons are made for rectangular wings over an aspect ratio, 
camber ratio and a range. 

Applications to rectangular wings of aspect ratio less than four 
are not attempted with the present method because only two spanwise 
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loading modes are employed. These two modes provide for only two 
spanwise rows of control points on a panel and are insufficient to 
insure tangent flow when the induced effects are large across the span. 
The use of additional spanwise modes from the same series (eq. (35a)) 
would allow for more spanwise rows of control points and hence provide 
better solutions. However, this is not attempted herein because of 
the substantial increase required in computer time. 

A study is also conducted using two-dimensional theory for a 
highly cambered section to: (l) provide limiting values for the three- 

dimensional results, and (2) establish relationships between the 
results at two different angles of attack which can be compared with 
those of the three-dimensional analysis . 

Variation of Results With- 


Box Size 

The integrations over the box of the approximated singular and 
nonsingular terms have yielded a result which is independent of the 
box size to 0 ( 5 °) . However, the size of the box can still have an 
effect on the answers by being either too large or too small. If the 
box is too large, it will include those portions of the wing which are 
outside the valid range of the functions approximating the singular 
part of the integrands. If, however, the box is too small the numerical 
integrations outside the box will suffer a loss in accuracy, unless 
substantially more integration stations are used, because of the close 
proximity of the edge of the box to the singularities. These two 
extremes are illustrated in Table I. 



TABLE I. 


Number of Integration Stations Required for Different Box Sizes 
and Effect of Box Size on Aerodynamic Characteristics for an 
A = 20, h/c = .00125 Rectangular Wing at a = 5° . 
Control Points at x/c = .3^5, .905* 


Box Size 

Integration Steps in 
Box Strip 

Aerodynamic Characteristics 


Chordwise 

Spanwise 

C L 

mm 

■SI 


■■ 

■ 

200 

51 

.5332 

.2521 

.0018 

-.0008 

HfiEi 

200 

51 

.1959 

.2559 

.0012 

.0031 


300 

151 

. 5121 

.2511 

.0011 

.0017 

EH 

500 

201 

.1450 

.3212 

.0031 

.0212 


Spanwise Control Point Locations 


Only two spanwise rows of control points are utilized along a 
semispan because only two spanwise loading modes are employed in the 
present analysis. A limited study of the effects on the aerodynamic 
characteristics of varying the location of these rows was made and the 
results, presented in Table II, show only slight variation with row 
location set. However, the 0.25, 0.75 set was chosen to be used 

TABLE II. 

Effect of Spanwise Location of Control Point Rows 
on Aerodynamic Characteristics of an 
A = 20, h/c = .00125 Rectangular 
Wing at a = 5° . Control 
Points at x/c =.345, .905 


2y/b 

C L 

cp 

C-p. . 

, 

c n . . 

D ,11 

o, .7070 

.1950 

.2562 

— 

.0031 

.25, .75 

.1961 

.2561 


.0035 

0 , .8165 

.5125 

.2557 

MM 

.0021 
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subsequently because it (l) resulted in the closest agreement between 
the two induced drag values and (2) because of symmetry represents a 
solution which constrains the flow at two additional control points on 
the left half of the wing as opposed to the total of three which 
result for a set with zero as one of the spanwise locations. 

The induced drag terms ^ and CD,ii 8X6 ca dd e d the f ar “ 
field and near-field results, respectively. determined by the 

method of reference 6, and from the combination of the 

distributed camber drag and the leading-edge suction as seen in appendix 
B. Exact numerical agreement between these two drag terms is taken to 
mean that the chord loading is correct and hence the best solution has 
been reached. In practice, exact agreement is seldom achieved, hence 
chordwise control point patterns which yield reasonably close agreement 
between the two drag terms are chosen subsequently (as above for the 
spanwise sets) as the ones to be used. 

Chordwise Control Point Locations 

For two control points, many placement patterns were investi- 
gated, as shown in Table III. From the table, it is evident that, for 
the first control point near x/c = 0.35, the variations of the aero- 
dynamic characteristics for the A = 5 , h/c = 0.00125 rectangular 
wing at a = 5° are reduced, especially those of the induced drag 
coefficient. By comparing the two induced drag terms in the table the 
closest agreement is seen to occur when the second control point is at 
x/c = 0.80 or 0.95. These locations are so near those prescribed by 
Multhopp in reference 6 for a two chordal loading solution, that is, 
x/c = 0.3^5 and 0.905 that his are used in the present method for 
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the two chordal loading solution. For higher camber ratios additional 
chordwise control points and loading functions are needed to effect a 
best solution as seen in the next section. 


TABLE III 

Effect of Locating Two Control Points on the Aerodynamic 
Characteristics of an A = 5, h/c = .00125 
Rectangular Wing at a = 5° 


Control Point Locations 

Aerodynamic Characteristics 

1 First 

Second 

C L 

X 

cp ^ 

wssm 

DH 

.i 

.9 

.3839 

.2398 

.0095 

.0060 

.15 

.80 

.3766 

.21+33 

.0092 

.0070 

.20 

.80 

.376U 

.21+1+6 

.0092 

.0071 

.20 

• 90 

.3803 

.21+62 

.0091+ 

.0073 

.25 

.75 

.3739 

.21+1+0 

.0090 

.0071 

.25 

.95 

.3823 

.21+79 

.0095 

.0075 

.35 

.65 

.3721+ 

.21+63 

.0090 

.0076 

.35 

.75 

.371+0 

.2I+7I+ 

.0090 

.0076 

.35 

.80 

.3759 

.21+85 

.0091 

.0078 

.35 

.85 

.3771+ 

.21+95 


.0078 

.35 

.90 

. 3781+ 

.2500 

.0093 

.0079 

.35 

.95 

. 3801 

.2511 

.0091+ 

.0081 


Number of Chordal Control Points 

The number of control points must increase with camber in order 
that the flow might be constrained to approximate better the normal 
velocity distribution or camber shape. If only two chordwise control 
points, hence loading functions are used, the aerodynamic characteristics 
which result would lose accuracy, especially those involving the 
pitching moment , because it depends more heavily on the third and higher 
modes than do the lift and leading-edge suction. 
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The effect of number and locations of control points on the 
aerodynamic characteristics is tabulated in Table IV for a typical wing. 
From this table it can be seen that there are several patterns which 
would appear to do equally as well from a best solution cheok. 


TABLE IV 


Effect of Number and Locations of Control Points on the 
Aerodynamic Characteristics of an A = 5, h/c = 0.031^ 
Rectangular Wing at a = 5° 


Control Point Locations 

Aerodynamic Characteristics 


Second 

Third 

Fourth 

Fifth 

C L 

x cp 

c n . 

D,i 

c D,ii 

. 3 U 5 

.905 




.8798 

.3761 

.0498 

.0077 

.05 

.35 

.95 



• 9936 

.3865 

.0633 

.0580 

.1 

.32 

.8 



- .0028 

13.6321 

.0000 

- .2656 

.1 

.32 

.9 



.5202 

• 5531 

.0176 

- 5.6964 

• IT 

.4 

• 95 



.5939 

.0773 

.0229 

- 2.0950 

.19 

.615 

.95 



.3877 

1.3917 

.0101 

- .9160 

.1*5 

.7 

.95 



.7685 

.1*995 

.0382 

.0504 

.52 

.78 

.98 



.81*82 

.1*038 

.0464 

.0361 

.05 

.35 

.65 

.95 


.8363 

.1*215 

.0451 

.0561 

.05 

.35 

.75 

.95 


.8331+ 

.1*223 

.0448 

.0561 

.1 

.2 

.8 

.9 


.5776 

.5331* 

.0218 

.0398 

.1 

.32 

.6 

.97 


- .3227 

. 817 I+ 

.0510 

- 1.5160 

.116 

.i*ii+ 

.75 

.97 


-3.11+21 

- . 3092 

.6342 

-68.5122 

.2 

.4 

.6 

.8 


.1356 

- 2.5608 

.0016 

- 2.2424 

.05 

.32 

• 55 

.75 

.95 

.8517 

.4296 

.0468 

.0583 

.08 

.29 

.57 

.83 

.98 

.8098 

.1+756 

.0424 

.0653 

J 

.32 

.5 

i 

• 7 

.9 

-3.8451 

.0277 

1.0037 

-71.4019 
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In order to try and determine if there was a preferential set of 
chordwise control points, a separate computer program was developed 
which solved equation ( b6 ) for the left-hand side using the previously 
determined { c nm K Hence, the variations from the tangency flow condition 
for points in between the control points could be examined and some 
indication of the best control point patterns obtained* From this 
study the pattern of x/c = 0.05, 0.35, 0.65, 0.95 was determined to be 
the best choice. It should be noted that this set may not be usable 
throughout the angle-of-attack range, for if one of these control points 
falls within the range 88.5° < 0 <91-5° a new point must be selected 
to avoid the numerical integration problems discussed in Chapter VII. 

In addition, this pattern may not be appropriate for wings of other 
camber or aspect ratio. 

Experimental Data Comparison 

Since the only reliable data available with which to compare 
was for flat rectangular wings, the present method, which is only valid 
for cambered wings, was studied numerically to determine what value of 
h/c would adequately represent a flat wing. A value of h/c = 0.000013 
was determined to be sufficiently small, as further small reductions in 
h/c didn f t change the numerical results. Two chordal loading modes 
with the associated control points already established were employed in 
obtaining solutions. The results of the computations, along with those 
of the modified Multhopp method (ref. 8), are presented in Table V and 
compared with those from the experiment on both a section and wing basis. 
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Further, all three lifting pressure distributions are given in figure 1 
at the plane of symmetry for an A = h rectangular wing at a = ^.35°. 

TABLE V 

Experimental and Theoretical Aerodynamic 
Characteristics for an A = k Flat 
Rectangular Wing 


Type 

2 y/b = 0 

: a = 4.35° 



C l 

X 

cp 

C L 

a 

X 

cp 

Experiment (ref. 2b) 

.323 

~=T 

CM 

.060 

.235 

Present Method (W=2) 
(2 semispan stations) 

.346 

.242 

.067 

.237 

Modified Multhopp Method (N=U) 
( 19 semispan stations) 

.331 

O 

CM 

.063 

.230 


Both the table and the figure show that the present method 
produces results which agree reasonably well with the experiment and 
the modified Multhopp method. Thus, having shown that the present 
method yields reliable results for flat wings , it is of interest to 
determine its efficacy for cambered wings. Because of the lack of 
reliable cambered wing data obtained for rectangular planforms 
meeting the camber constraint of the present method, as mentioned 
in the Introduction, cambered wing correlations are made with 
results obtained with the exact two-dimensional method given in 
appendix C. 














Figure I.- Lifting pressure coefficient distribution for A = 4 flat rectangular wing 
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Theoretical Comparisons 

Two- and Three-Dimensional Lifting Pressures 

The lifting pressures for the two-dimensional transformed 
solution (appendix C) are compared with those of the present method 
and modified Multhopp method at the plane of symmetry of an aspect 
ratio 20 rectangular wing, with all having h/c = 0.00125 and at 
a = 5°. The results are graphed in figure 2(a) and show that the 
present method for N = 2 agrees better with the two-dimensional 
than N = 3 for the present method or N = 2 for the modified 
Multhopp method. A similar graph is presented in figure 2(b) for 
h/c = 0.0314 and a = 5° with N = h for the present method and 
N » h for the modified Multhopp method. A comparison of the three 
sets of data shows that the present method predicts better the two- 
dimensional c x value, whereas the modified Multhopp predicts better 
the x C p value and lifting pressure distribution. The poor 
quality of the lifting pressure distribution predicted by the present 
method is serious and is discussed later. 

Two-Dimensional Cambered Wings 

The prediction of the aerodynamic characteristics for two- 
dimensional circular-arc-cambered wings can be accomplished in either 
of two ways. The first is in the classic manner of transforming a 
circle into the airfoil shape and arriving at the moment and forces by 
Blasius* theorems as given in appendix C (referred to herein as the 




k3 



(b) h/c = 0.031 1 * 
Figure 2. - Concluded. 


Xc P 

.355 

.363 

.315 


kb 


transform method.) The second is based on thin airfoil theory (small 
a and h/c approximations) and employs the Birnbaum series to 
represent the chord loading and from which the moment and forces are 
determined (Reference 25 has examples of the procedures.). When the 
results of these two methods are compared in Table VI for h/c = 0.50 
and a values of 0° and 5° , remarkable agreement is noted for 
both the x C p and values. In order to understand this 

agreement, the chord loading predicted by each method at a = 5° is 
graphed and appears in figure 3. An examination of this figure shows 
noticeable disagreement in lifting pressures along the chord; however, 
the integrated effect of these differences is seen from Table VI to be 
small because of the compensating behavior of the disagreements. 

TABLE VI 

Circular -Arc Lift and Center of Pressure 


— 

i 

Method 

c^ Equation 

x Equation 

cp 



Transform 

27r[sina + 2(h/c)cosa] 

sina cosa 

2 - sina - 

C0Sa (f ") 



C^TT 


Thin airfoil 

_ /air _ h\ 

2 niBo + 2 o) 

1 . TT(h/c) 
F + c 

1 











Transform i 
Thin airfoil 



Figure 3. - Two-dimensional lifting pressure coefficient distribution 

for h/c =0.5 at a = 5°. 
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TABLE VI continued 




h/c = 

.5 


Method 

a 

o 

o 

II 

a = 5 

O 


c 

X 

c 

X 


1 

cp 

i 

cp 

Transform 

6.2832 

.50 

6 . 8069 

.1+800 

Thin airfoil 

6.2832 

.50 

6.8315 

.1+799 


Note, in particular, that the center-of-pressure predictions hy both 
procedures yield the same results at both angles of attack; whereas, the 
thin airfoil theory predicts a higher lift than the transform method 
only at a = 5°. This relationship for the lift at a > 0° between 
the solution on the surface (transform) and along the chord (thin 
airfoil) will be useful later. 

Over an Aspect Ratio Range 

The effects of aspect ratio on the center of pressure can be 
seen in figure 4 for both cambered and flat wings. It is interest- 
ing to note that for decreasing aspect ratio the predicted flat wing 
centers of pressure move forward as expected; whereas, for cambered 
wings at a = 0° the centers of pressure progress rearward when 
computed by all three theoretical methods. A comparison of the curves 
shows that the present method and the modified Multhopp cambered-wing- 
x C p results agree well and tend toward the two-dimensional limit. They 
both also show a more rearward location of at all aspect ratios 
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X Cp 



— CAMBERED WINGS 

x Two-dimensional 

Present method 

Modified Multhopp method (ref. 8) 

Scholz (re f. 23) 

/ FLAT WINGS 



.5 X 


4 


.3 


I/A 


Figure 4. - Effect of camber on center of pressure. Cambered wing 

results are at a = 0°. 



U8 


than those predicted by Scholz. For the flat wings the present method 
results tend toward the two-dimensional value of x C p slightly faster 
than the other two methods with increasing aspect ratio. At a = 0° , 
the effect of camber cannot be seen at a given aspect ratio because the 
different camber curves tend to collapse upon one another. Hence, these 
effects are discussed in the next section at a given aspect ratio and 
positive angle of attack. 



Sketch 9 - Loading composition for finite-aspect ratio 
flat and cambered wings 

The above center-of-pressure trends can best be understood by 
examining the preceding sketches. The effects of finite aspect ratio are 
seen to induce both a negative cot e/2 load and a negative sin e load 
with these induced loads becoming larger as the aspect ratio decreases. 
For flat wings, the induced cot e/2 load does not change the load 
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center but does decrease the lift at a given a ; however, not only does 
the sin e load decrease lift, it also imposes a nose-up moment about 
the two-dimensional lift center as well as the leading edge. These 
effects lead to a forward movement of the center of pressure with de- 
creasing aspect ratio. For symmetrically cambered wings (in particular, 
circular-arc wings at a = 0°), the induced sin e load reduces the lift 
at a given a without distrubing the center of pressure; whereas the 
induced cot e/2 load decreases the lift and produces a nose-up moment 
about the leading edge or a nose-down moment about the two-dimensional 
lift center. These two effects cause the center of pressure on cam- 
bered wings to move increasingly rearward with decreasing aspect ratio. 

The trend of center of pressure with aspect ratio is well 
known for flat wings; however, the corresponding trend for cambered 
wings is not as well known, even though documented by Scholz in 
reference 23 • 

Over a Camber Ratio Range 

The effects of camber ratio on x C p and are seen in 

figures 5(a) and 5(b), respectively, for an aspect ratio 5 rectangular 
wing at a = 5° • As expected, the predicted results show an increase 
in and rearward movement of x C p with increasing h/c . Further- 

more, comparisons made between the modified Multhopp method and the 
present method indicate that the centers of pressure agree reasonably well 
up to h/c of .013; whereas, the predicted by the present method 

exceed those of the modified Multhopp method for all camber and increases 


at a faster rate. 



x C p 

.3 


.2 

Modified Multhopp method (N = 4) 

Present method (N = 4) 



Figure 5. - Effect of camber on some aerodynamic characteristics of an 
A = 5 rectangular wing at oc = 5°. 
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Over an Angle-of-Attack Range 

The variation of with a at an aspect ratio of 5 and 

camber ratio of 0.031^ are computed for both the modified Multhopp and 
present methods, and the results presented in figure 6. The curves 
show that the present method predicts a higher value of at all 

angles of attack and yields a higher lift-curve slope. 

Discussion 

From the preceding studies of the present method and the 
subsequent applications, much has been learned: 

(l) Number and Locations of Control Points 

Two control points (loading functions) are sufficient to 
describe adequately the wing having only small amounts of camber 
(h/c - 0.00125), and the results with different locations indicate a 
low sensitivity to position. This happens because the first two 
lifting pressure functions are adequate to describe the pressure 
loading and also satisfy the tangent flow boundary condition all along 
the arc. For higher camber ratios, more than two control points are 
required to constrain the flow sufficiently to meet the solution 
criteria. The number and location of these control points depend 
upon wing aspect ratio, and angle of attack because of the few control 
points constraining the flow. Numerical studies indicate that two 
additional terms in the Birnbaum pressure loading series are still not 
adequate to describe the loading distribution for wings with cambers in 
excess of 1.25$; therefore, a larger number of terms are required. 
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This is not a very practical solution since it necessitates an excessive 
amount of computer time. For example, with two spanwise stations and four 
chordwise loading functions , computer time on the CDC 6600 of 15 minutes 
is required, and for six chordwise modes, time in excess of 30 minutes 
would be needed. Ideally then one would utilize a better set of modal 
functions. The problem is the determination of this better set. One 
procedure for doing this would be to obtain them from the local pressure 
equation given in appendix C for the exact two-dimensional solution. 
Another would be to obtain them from the two-dimensional integral 
equation relating the downwash ratio or local slope to the pressure 
distribution over Ax . Still a third procedure would be to express the 
series of loading functions in terms of distance along the arc having the 
appropriate singularity as the first term with the other terms having 
the first term multiplied by distance along the arc to a power. 

The first procedure would be very difficult, as an examination 
of equations C-23 and C-2k will easily show. The second procedure 
leads to another series solution which is as yet undetermined. It is 
likely, however, that many terms would be required. The third procedure 
has been programmed , with only slight improvements obtained. 

(2) Limiting Values 

The present solution tends to the correct two-dimensional value 
of x for both the flat wing and at a = 0° for the cambered wing. 

The present method also predicts a C-^ which tends to the two-dimen- 
sional value for h/c - 0.031^ and a = 5° • It is interesting to note 
that, whereas the two-dimensional solution along the chord (thin airfoil) 



predicts a higher c^ than the solution along the surface (transform), 

the reverse is true for the two three-dimensional solutions at both 

A = 20 and A = 5 • 

(3) Moderate Aspect Ratio Comparison 

The present method appears to predict reasonably well the 

aerodynamic characteristics of wings with small camber (h/c - 0.0125). 

This is demonstrated by comparison with the experiment and the modified 

Multhopp method. The fact that at h/c = 0.031^ the present method yields 

a higher C T at any a , a higher Cr , and a more rearward x 

ij a v 

than the modified Multhopp method can be attributed to both the modal 
functions employed in the present method and the limitations inherent 
in the planar approach. It should be noted that even with the modal 
functions employed satisfying the solution criteria, poor results can be 


obtained because of the limited number of terms used. 



IX. CONCLUSIONS 


The development of a nonplanar lifting surface method having a 
continuous distribution of singularities and satisfying the tangent flow 
boundary condition on the mean camber surface is given in this 
dissertation. The method predicts some incompressible longitudinal 
aerodynamic coefficients of rectangular wings which have circular-arc 
camber. After some preliminary sensitivity studies are conducted with 
this method, applications are made over an aspect- and camber-ratio and 
angle-of-attack range and the results compared herein with other 
theoretical methods and flat-wing experimental data. From these studies 
and comparisons, several conclusions emerge. They are: (l) The present 

method is able to predict with good accuracy some longitudinal aerodynamic 
characteristics for camber ratios not greater than 0.00125 at both high 
and moderate aspect ratio by employing only two control points (with 
little sensitivity of results to control point locations) and without 
the singular strip encountered in the flat-wing theoretical solutions . 

(2) For higher camber ratios, more control points are required, but even 
with their locations yielding results which give reasonable agreement 
between the far-field and near-field induced drag values (the solution 
criteria) , the lifting pressure distributions can be poor. The poor 
distributions and high sensitivity of results to control point location 
indicate that the pressure mode functions chosen are not the most 
appropriate for wihgs having cambers in excess of 1.25$ • Other lifting 
pressure functions were investigated without obtaining any significant 
improvement in the results. (3) The control point locations are 
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dependent upon aspect ratio (because of the few control points used) , 
camber ratio, and angle of attack, because of the inadequacy of the modal 
functions chosen. (4) The present method is also able to predict 
reasonably well the lift coefficient and center of pressure for moderate 
aspect ratio (aspect ratio 5 illustrated) at an angle of attack of 5° 
up to a camber ratio of about 1 . 25 # • (5) The present method’s 

predictions at a moderate aspect ratio of lift coefficient at any angle 
of attack and lift-curve slope, as well as center of pressure, are larger 
than those of the modified Multhopp method for all camber ratios. The 
over-prediction of the lift coefficient by the present method when 
compared with the modified Multhopp method is the reverse of that found 
in comparison of the two-dimensional surface (transformed) and chordal 
(thin airfoil) solutions. This indicates that the modified Multhopp 
results are low for certain configurations and that for others the 
lifting pressure functions of the present method are not appropriate. 



X. RECOMMENDATIONS 


The work reported herein could serve as a model for lifting 
surface solutions of other basic camber shapes such as elliptic, 
parabolic, or those represented by some other general function. 
However, it should be kept in mind that each of these other camber 
shapes has a different set of appropriate chordal loading functions 
which must be determined in order to yield acceptable results in a 
reasonable amount of computer time. Another way to save computer 
time would be to perform the general spanwise integration analytically 
first, leaving only the chordwise integration to be done numerically. 

Extensions to the present work could be undertaken to include 
both the effects of sweepback, without introducing reversed twist, 
and variable dihedral. 
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XIII. APPENDIX A 


Certain Integrals Over The Box 


Some of the integrals presented in equation (4l) for integration 


over the box contain in the integrand the term 


’COS 2 0+((jl)/R) ‘ 


solution for these definite integrals can be analytically formulated; 
however, obtaining correct results from the numerical evaluation of 
these formulations may be difficult. The difficulty is evident for 0 
near tt/2, because there the numerical evaluation of the general formu- 
lation does not tend smoothly to the same numerical result as would 
be obtained by a direct integration of the original integrals for 
0 = tt/2. The numerical difficulty arises because of the general 
occurrence of (l/cos0) as an overall multiplier, as can be seen in 
the following example. 

Example : 

Consider the surface integration of 


f° f R q 2 d((jq/R)da 

J- M [c 2 cos 2 0 + (w/R) 2 ] 2 

— o ^ 


(A-l) 


-a 

R 


Integrating expression (A-l) with respect to w/R and imposing the limits 


leads to 


R 

COS ^0 


(u)/R)da 

a 2 cos 2 0 + (w/R) 


[ ( U ) 

I tan~ 1 \Racos0 / 

J - o 

-a J 


(A- 2.) 


The first integral of expression ( A-2) can be evaluated readily; how- 
ever, a change of variable is needed in the second. These steps 
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lead to 


RacosS 


R 

2tan _1 ( Rac _°s0 \ 

a) 

J. f Ftt/ 2 - tan _1 (x)l 

(A-3a ) 

cos 0 

\ W / 

J x 


where 


-R0COS0 

m 

x = Racos0/(jj 

( A-3b ) 


The first term in the remaining integral integrates to zero because 
it is odd over the integration range. The second term has to be 
integrated term by term after expanding tan -1 (x) into 


tan _1 (x) = x - — 2— + — 2— 2— + ... 

3 5 7 

Hence, upon performing the above integrations and letting 

Rj = w 

leads to 

2tan -1 ( cos 6) - 2 (cos0 - cos3e . cos5e cos7e + . . . Y| 

cos 3 0 I l 9 + 25 49 J 


(A-4) 


(A-5) 


(A-6) 


This expression could be considered as the general formulation 
of the evaluated definite integral appearing in expression (A-l.) . Note 
that, in expression (A-6)» values of 0 near tt/2 produce large results 
for the (l/cos 3 0) multiplier and small results for the terms inside 
the braces. Now, upon performing the indicated multiplications, 
there result some terms which are still small and tend to zero as 0 
goes to tt/2; whereas, other products, 2Rtan ^(cos0) -2Rcos0 

cosre cos 3 0 5 

produce numerical indeterminancies or, at best, infinities of opposite 
sign. In either case, the answer or finite part of the surface 
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integral is not determinable. 

This problem can be circumvented by expanding in expression 
(A -6) the term tan*" 1 (cos0) according to equation ( A— U) and combining 
like terms inside the braces. After simplifying, the result is 

H ^ 1 cos 2 e - cos 4 e + • • • J (A-T ) 

Now, for 0 = tt/2, the expression, instead of producing indeterminancies, 
gives a value of - \ R. This result agrees with the direct integration 
of expression (A-l) for 0 = tt/2. 

Another integral which must he integrated in the same manner 
as the preceding is 


R 



(u/R) 

(w/R) 


d(qj/R)aq 

a 2 cos 2 0 + (w/r) 2 


(a-8) 


and results in 

-^R(l - ^cos 2 0 + “cos 4 0- . . ) 


(A-9) 



XIV. APPENDIX B 


Determination of Aerodynamic Coefficients 

The determination of some section and wing aerodynamic 
coefficients for wings with circular-arc camber is accomplished by using 
the c^ set computed from equation (^6) and the equations given in the 
following paragraphs . 

Section Coefficients 

The computational equations for the section lift (including the 
contribution due to leading-edge suction) pitching moment, distributed 
camber drag and leading-edge suction coefficients are given below. 


Section lift coefficient: 


c (y . ) = — 
l 7T 


~om 


2tt sin 2 9 


V !£• rh n (e)»irf!i 

tkX&o^'o * 

cos6 ° ,| p Tv/n 

o sin(0 o + a)J\b/2/ V 


inede 

al 


(B-l) 


Section pitching-moment coefficient : 
2 N-l 


2 N-l rn- 

v ) = ? f 2 2 Cnm j ' 

LE m=o n=o ° 

even 


h (e)sin( ■O - 0 )sinedeAr, 
n o 


sin($ + a) 


\ h/2 / 


i- 


f y. 


(B-2 ) 


Section distributed camber drag coefficient : 


^ r 7r h n (e)cosT3sin; 

^ tt c nm| sin( t 9 + a 

m=o n=o •'o 


7r h^(e)cosT3sin£de /y . 

PVV2J 


1 - D4- 


y j 

b/2. 


(B-3) 


even 
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Section leading-edge suction coefficient : 

2 


LE 


/ 4 i 4.i \ m L /y.i V 

^5^ 2tt sin 2 0 o sin(0 o + a) v 5 / / * \ ^ / 

even 


(B-U) 


where 


•3 = cos" 1 ^ 1 ~^ c . P s£ ^ £ + cos(0 q + a)] -a (B-5) 


Wing Coefficients 

The equations used to determine the wing lift, pitching moment, 
distributed camber drag, leading-edge suction and near -field induced-drag 
coefficients were determined by Multhopp quadrature and axe given as 
follows for symmetrical span loadings . 


Wing lift coefficient: 
Jvl 


C— = 

L c 


: av (/ + iy (S C i (y j } c( yj ) sin ^ C ( y M) 


(B-6) 


Wing pitching-moment coefficient : 

J-l 
2 


TT 


m LE c c ,( J + 1] 
av ref 


n5 '“ls (yj) cHv3) Sin * 3 + * \ E ( y iti) 


(B-7) 


Wing distributed camber drag coefficient : 




C D = c C 
c av 


tVtt( 2 c(yj) sin^ + | c d c ^jrtj 


(B-8) 
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Wing leading-edge suction coefficient : 

j _ 1 


LE 


1) (^' 1 CS LE <y f ) ° <y ^ Sl "^ + 2 'fe) 


2 

(B-9) 


Wing near-field induced-drag coefficient : 


C D.. = c d - C Stt? sin0 o 

li C JL»Xli 


(B-10) 


where 




(B-ll) 


Wing far-field induced-drag coefficients : 


: D. = A f 5iT “i dr > 
1 •'-1 


(B-12) 



XV. APPENDIX C 


Two-Dimensional Circular -Arc Airfoil 
Aerodynamic Characteristics 


It is well known (see refs. 26 to 29) that a circle in the complex 
C-plane can be mapped into a circular-arc airfoil in the complex z-plane 
by selecting the origin of the circle to lie at ih/2 and then applying 
the Kutta-Joukowski transform 


z 



(C-l ) 


to points on the perimeter of the circle (see sketches 10 and ll). The 



circle and circular-arc airfoil. 
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Sketch 11*- Additional geometrical relationships between 
generating circle and circular- arc airfoil. 


points from counterclockwise to will lie on the upper surface 

and those from counterclockwise to will lie on the lower 

surface. The local velocity, and hence pressure, can be determined over 
the airfoil by making use of the complex potential (w) of the circle of 
an angle of attack and the Kutta-Joukowski transformation as follows : 


dw _ dw dC 

dz dC dz 


(C-2) 


where 


w = V£e + 


, r 2 -ia 
ia Va e 


C - s 


+ 2aiVsin3 ln(t; - s) 


(C-3) 




TO 


s location of the center of the circle in 5-plane 


3 the angle , where is the point on the 

circle that transforms into the trailing edge of the 
airfoil (point H ) . 


The angle 3 must he such that the velocity at the point H remains 
finite. To solve for the required 3 relationship, the procedure 
given in reference 2 6 will he followed. Begin with 


dw , r ia 


Va 2 e~ ia 2aiVsinB 
(5 - s) 2 5 - s 


(C-Ua) 


or re-express it as 


ff=4 ia 



ae- i(a + P) \ 
5 - s / 


and the differential of equation (C-l) 


dz 



(C-Uh) 


(C-5) 


Upon substituting equations (C-Uh) and (C-5) into (C-2), it follows that 

a -i(a + 3)' 

5 - s / \ r z _ ^ 


dw 

dz 


- 4 la ♦ i^-X 1 - 


ae 


^ ) 


(c-6) 


The point on the circle can be seen in sketch 11 (ref. 27) and 

corresponds to £ = -k . In order to keep the velocity from taking on 
an infinite value at the transformed point H , the point must be 

made a stagnation point (dw/d£ = 0) so that dw/dz -will assume 
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the indeterminant form of 0/0 , which can be shown to yield a finite 
velocity at the trailing edge. 

Next, the angle 3 must be related to other known quantities. 

This is accomplished by examining the geometry of the generating circle 
in sketch'll. From it the distance along £ from to the 

horizontal projection of the circle center is determined to be 

k + = acos(3 - a) + 0 = acos3 (C-7a) 

Furthermore, the vertical distance along ir| is found to be 

s = asin(3 - a) = asin(3 (C-7b) 

n 

These can be written as 

k + is^ = a[cos(3 - a) + isin(3 - a)] = ae^^ (C-8a) 

or 

ae"^ = (k + s) e^ a (C-8b) 


By making use of equation (C-8b) in the first group of terms in equation 
(Ct 6), the following result can be obtained: 


dw 

dz 


■ 4 io + - 


ae 


-i(a + p) 


? - s 


) (C - + k) 


(C-9) 


which after some simplification can also be given as 


dw ir ia 
— = Ve 
dz 


(C - s)- ae 


-i(a + 3) 


U - s) 2 (? - k) 


(C-10) 
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With the addition and subtraction of k to the bracketed term, and upon 
using the result of equation (C-8b), equation (C-10) can be recast as 


dw ia/_ , , i(3 - a) 

— = Ve I 5 + k - ae - ae 


-i(a + 3)j 


(C - s) 2 (S - k) 


(C-ll) 


By expanding the exponent terms, equation (C-ll) becomes 


|| = Ve ia ( ? + k - 2ae- ia coss) 


(; - s) 2 (C - k) 


(C-12 ) 


An examination of the velocity at the point H (or where £ = -k) yields 


U H “ 1V H ~ 


dw , r 1 Q!/ a --LU n \ 

- — = Ve (0 - e cos3 ) 
dz 


ae i2(B - o) 


(0-13) 


or 


U H * 1V H = " 


ae 


Vkcos3 
i(23 - 2a) 


<c-H0 


In order to obtain the conjugate of u - iv the angle 23 - 2a can 

n n 

be rotated tt radians. That is 


U H + iV H = 


Vkcos3 i (23 - 2a + tt) 

e 

a 


(C-15) 


Hence, if the angle that a tangent to the cusp at the trailing edge makes 
with the x-axis is 23 - 2a, finite velocities are obtained. 

From sketch 11 ? it can be seen that 


a 2 = k 2 sec 2 3 


(C-l6) 
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and 


k = c/4 


(C-17) 


When these are used in combination with equation (C-8b), the following 
can be written: 


!««» •**•(% *£)•*“ 


(C-l8a) 


Therefore, 


/f + it 

e« 

sec 3y 


(C-l8b) 


Taking the natural log of both sides leads to 


i2h 


i3 = ln^cosp + cosp/ + ia 


3p) 


(C-19a) 


or 


ln(cos 2 P + — cos 2 P ) 

ip = — - 1 + itan 


■a 


+ ia (C-19b) 


which yields 


- tan 1 + a = 3 + a 


(C-20 ) 


In order to determine the pressures on the airfoil, it is more 
convenient to specify locations on the circle and then determine the 
corresponding point on the airfoil rather than conversely. Hence, by 


defining 



(C- 21 ) 


Tfc 




ae 


iY 


+ 


ih 


and by substituting this along with equations (C-l6) , (C-17) and (C-20) 
into equation ( C-12) } the following equation can be written: 

dw „ iaf iy . ih . c 2k -ia /7 r . N \ 

55 = Ve V ae + ~ IT " 555e e cos(|3 * a) )' 


Pllgl* 


(C-22) 


After all the multiplications and simplifications have been carried out, 
the velocities are determined from 

dw Vcosg||cosaP^ ~ sinaPgj + i jcosaPg + sinaP-J| 

(I)[i - cos (y - e)] 

(C-23) 


U - IV 


where 


dz 


= -acos 


(y + 3) + |a - o) r cos( 3 + hsing + a)^sinp| 


+ cos V[% + I^ cos0 + a 

ho)^ 


hai^sinP hco^. 


a 

+ siny|o3 T + - 7 *- sin(3 + h + sin(3 - — cos (3 


siny |o 


I a 


B *2 w ph 
. 2.1 h . R 


* cos ^[- 55 + it 2- cosB - vr - 55 

+ sin 2 y 


hu)^ h 00 


cos 3 j 

^-*j 

5 ) 


7 — sing 


fh 2 u 

ho: 

h 2 w 1 

1 

cos(3 + + 

— sing 1 

Ilia 2 

a Ua 2 J 


rh 2 u. 


+ cos3y -I 


R 


■ 4a 2 


+ sin3y 


1 4a 2 J 


(C-23a) 
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and 


- -asin(y + p) + [-w^cosp - hcosp - w^sinp] 

hui^cosP 1 


+ cosylWj + — sinp + h + sinp - 


] 


-siny 


[ , 2 hw_ hu T 

Y ’ SI cos5 * — 

I 


K _ hw h V 

+ cos2y| cosp + + sinp 

Ua ! a lta 2 

. , [ h 2 h “R h Ul h 2 ^ 1 

-sin2Y ■&+— =°se - — - — 

2 ,. 


r h 2 w l [ h 2 w ] 

cos3y - - sin3y| 1 

Ua 2 J 1 Ua 2 J 


(C-23b) 


“R 


co T = 


— seep cosa cos(P + a) + -jj- 


h c 

— + — seep sina cos(P + a) 


(C-23c ) 
(C-23d) 


Hence 



(C-2Ua) 


(C-2Ub) 


The location on the airfoil which corresponds to the points on the circle 
at which the velocities sire determined can he computed by substituting 
equation (C-21) into equation (C^l) and simplifying the resulting 
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expression. The result is 


-*[(' 
i 


sec3 cosy + 


cos3 cosa 


,'„2o) 


+ i(sec3 siny + 


1 + 2siny sin3 + sin 2 3 

tan3 - - c - ^ (sin T t_sing) \] 

1 + 2siny sin3 + sin 2 3'-* 


(C-25) 


These equations have been programmed and the lifting pressure results 
are compared with the three-dimensional solution at the mid-span of a 
rectangular wing with a very large aspect ratio in the Results and 
Discussion chapter of the dissertation. 

The forces and pitching moment which are developed over the circular 
arc airfoil can be determined by direct integration of the pressures and 
the moment of the pressure over the airfoil or by employing the analogous 
integrations which result from Blasius' theorems. The Blasius' Theorem 
for forces (ref. 26) can be expressed as 


X - 



dz 


(C-26) 


or. 


X - iY = 



(C-27a) 


By using equations (C-Ub) and (C-5), equation (C-27a) can be written as 

x . s)e la * 

.[(C-s)-ae- i(c, + P) ] 2 


<C - - k 2 ) 


(C-27b) 



After expanding all of the terms and employing the Residue Theorem, the 


result is 


Hence 


X - iY = -i87Tq asin3e ia 


and 


X = q 87ra sing sina 


Y = q 877 a sin3 cosa 

With the substitutions for a , equation ( C-l6 ) , and 3 
the coefficients of forces in the x and y directions 

c x = 27T(sin3 cosa + cos3 sina) sec3 sina 


or 


and 


c 


x 


= 2tt| 


sin 2 a 


+ sina cosa 



or 


Since 


c = 27r(sin3 cosa + cos3 sina) sec3 cosa 


c 

y 


27r|sina cosa + 


cos 2 a 



c = c sina + c ir cosa 

\ A J 

The resulting equation obtained by substituting in for 
equations (C-30b) and (C-30d), respectively, is 


( C-28) 

(C-29a) 

(C-29b) 

, equation (C-20), 
can be written as 

(C-30a) 

(C-30b) 

( C-30c ) 
(C-30d) 
(C-31) 

; x and Cy the 
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Also , 


c d = c x cosa - Cy sina 

and for the same substitutions as before leads to 


c d = ° 


(C-32b) 


(C-32c) 


The Blasius Theorem for pitching moment from reference 2 6 about the 
center of the generating circle is simply 

“[-Hf-’cKifH (c - 33a) 

After substituting in it the results of equation (C-l) and sketch 11, 
this equation becomes 

“<= - R [- b /(« ♦ r - r)(;|) a (S) «l , (c - 33b) 


and with the inclusion of the results of equations (C-4b) and (C-5)j 
equation (C-33b) can he written as 


M C = R 


[- - f\ l (? - s)eia * ael6 ] 2 * 

\u - .) -ae- 1 '" - »>]’ 


(C - s) 4 (? 2 - k 2 ) 


(C-33c ) 


Upon expanding these terms, combining, simplifying and integrating each 
set of terms by the Residue Theorem, the following equation results: 

c 2 

M c = q ttjj— sin2a 


(C-3Ua) 
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or 


c M = i sin2a 

c 


(C-34-b) 


After expressing a in a manner similar to reference 28, that is, 


a = 0 + (a - 3) 


(C-35) 


making use of the identity that 

sin2a = sin2[p + (a - 0)] = -sin20 + 2cos(a - p)sin(a + 0) ( C— 36 ) 

and re-expressing (C-32a) in terms of 0 , 

= 2irse«0[sin(a + 0)] (C-37a) 


or 


_ c COS0 
sin(a + 0) = — 1 4 


2tt 


(C-3Tb) 


equation (C-3^h) can he rewritten as 


TVL 


* f[- si 


sin20 + 


2cos(a - ff)c cos0 
2ir 


( C~38a) 


or 


7r - c i - 

C M = " ]j" s i n2 & + ]“ cos3 cos(a - 3 ) 


(C-38b) 


The first term is split off and defined by 


7T 


C M ■ “ IT Sin2 P t 


(C-39) 
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which is the pitching moment about the focus and is independent of a . 
The pitching-moment coefficient about the leading edge is obtained by- 
transferring the moment about the center of the circle to the leading 
edge. Hence 




sin2a - 2 tt 


sina + cosa 



(C-UOa) 


or 



[ sinacosa 
2 


- sina - cosa 



(C-AOb) 


As a consequence of equation (C-40b) and equation (C-32a), 
■written as 


x 

cp 


can be 


x 


cp 


sjjaacosa _ 

2 


sina 


-2 si 


sina + cosa 


- cosa 

w 



(c-ia) 


These equations have also been programmed and results are compared with 
the three-dimensional solution in its limit in the Results and 


Discussion chapter of the dissertation. 



